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r-| ! Abstract. 

We have extended the hyperspherical variables method to the analytical calculation 
of the angular integral of the box graph. We discuss the applications of our results to 
the analytical calculation of the QED contribution to the electron g-2 of the set of 
three-loop triple-cross vertex graphs. 
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1. New planned high precision experimental determinations of the electron g-2 demand 
correspondent higher precision calculations of QED theoretical contributions. One-loop, 
two-loop and almost all the three-loop QED contributions to the electron g-2 are known 
in analytical form. The analytical calculations of the known three-loop contributions 
relied on two approaches: the hyperspherical variables method [1] [2] and the dispersion 
relations method [2] [3], the choice being dependent on the topology of the graphs. A 
hybrid dispersive and hyperspherical approach turned out to be more convenient for 
some troublesome sets of graphs [4] [5] . 

There are two sets of three-loop graphs whose contributions are still not known 
in analytical form, but only in numerical form. Analytical calculation of one of them, 
the set of graph called "corner-ladder" , is in progress [6] . On the contrary, none of the 
above mentioned analytical approaches seems to be suitable for the analytical calculation 
of the last remaining set of (not planar) graphs called "triple-cross" (see Fig.l). The 
application of the dispersion relations method to the whole graph seems difficult because 
of the high number of cut graphs and the complexity of the calculations involved. On the 
other hand, the hyperspherical variables method (in the formulation of ref.[l]) cannot 
be applied to box graphs or not planar vertex graphs at all. No useful way of combining 
both methods has been found so far. 

In this paper we present an extension of the hyperspherical variables method to 
the box graph and we discuss the application of our results to the calculation of a 
contribution of a typical not planar "triple-cross" graph. 

2. Let us summarize the main steps of the analytical calculation of a Feynman integral 
by means of the hyperspherical variables approach [1] : 

1) The Feynman integral is Wick-rotated to an Euclidean integral; so doing, the 
external momenta become all spacelike, and the result of the integration must be even- 
tually continued to the timelike region. 

2) Four-dimensional hyperspherical polar coordinates are introduced for each loop 
momentum. 

3) The angular integrations are analytically performed. 

4) The remaining radial integrations (one for each loop) are performed analytically. 

Obviously the points 1) and 2) do not depend on the particular kind of graph con- 
sidered. The point 3) is more critical; an elegant and fast method for performing the 
angular integrations is presented in ref.[l]: the denominators are expanded in Gegen- 
bauer polynomials (related to four-dimensional hyperspherical harmonics), the angular 
integrals are performed using orthogonality relations and the infinite series so obtained 
is summed up analytically. The feasibility of the subsequent radial integration is re- 
lated to the analytical form of the results of the angular integrals. As an example, let 
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us consider the simple 4-dimensiorial angular integral of a scalar bubble-loop graph 



*k (i) 

(q-k) 2 + m 2 ' 1 ' 

where q and k are vectors in a 4-dimensional Euclidean space. The denominator of 
eq.(l) is expanded in Gegenbauer polynomials C^(q ■ k): 
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(q-k) 2 + m 2 \q\\k\l-2Zq-k + Z 2 
where 



1 oo 



_ q 2 + k 2 + m 2 - R(q 2 , k 2 , -m 2 ) 

~ 2\q\\k\ 



R(x,y,z) = a/x 2 + y 2 + z 2 — 2xy — 2xz — 2yz ; 
the angular integrations are performed by using the relation 



J dQ q C 1 n (q-k) = 27T 2 S 0n , 



so that the summation of the infinite series (2) becomes trivial and 

(g — k) 2 + m 2 \l\\k\ 

Successful application of the above described method to more complex angular inte- 
grals depends strongly on the momentum routing and of topology of the correspondent 
graphs; in fact, these conditions must be satisfied [1]: 

i) The momentum of each line of the graph must be a linear combination of only 
two distinct momenta, loop or external; for every planar vertex graph one can find a 
momentum routing which satisfies such a condition, while in the case of the not planar 
"triple-cross" vertex graphs this is not possible. 

ii) At least one vector must not appear in more than two Gegenbauer polynomials, 
as tables of integrals of product of three Gegenbauer polynomials containing the same 
vector are not known; the box graph and all graphs made up only of box subgraphs do 
not satisfy this condition. 

Let us now consider the angular integral containing two denominators 
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d% (4) 



[(q - h) 2 + m 2 ][(q ~ k 2 ) 2 + m 2 } 
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which is the hyper spherical angular integral of a triangle loop; this integral can be 
worked out in the same way of previous integral (1); in fact, applying the expansion (2) 
to both denominators, one works out the double series 



/ 



dtt q 1 

x 



[(<? - ^i) 2 + m i][{Q - k 2 ) 2 + m%] <? 2 |/ci||/c2| 

oo oo 

x E E Z r +1Z 2 +1 Cl (q ■ h) Cl (q ■ k 2 ) , (5) 



m=0 n=0 

q 2 + k 2 + mf - R(q 2 , fc?, -mf) 



Zi = 



2\q\\ki\ 

which, by using the orthogonality properties of the hyperspherical harmonics 

J C l m {q ■ h) CHq ■ k) dn q = 2tt 2 ^fjC^Cki ■ h) , (6) 

is reduced to one single series whose sum contains two logarithms of complex arguments 
containing Z\ and Z 2 ; we do not show this result here. 

On the contrary, we underline a simple but important fact: the angular integrals 
(1) and (4) can be calculated directly integrating over the three hyperspherical angles, 
avoiding the use of the expansion (2), and writing 

dQ q = sin 2 M#i / sin6 2 d9 2 d$ 3 ; (7) 
Jo Jo Jo 

the calculations become rather cumbersome, but, as it will be shown later, they are still 
feasible when the expansion (2) cannot be applied. 

Let us now calculate the integral (4) in analytical form using the decomposition 
(7). We combine the denominators by using one Feynman parameter, the simple angular 
integration is done by using the 1-denominator result (3), and the subsequent integration 
over the Feynman parameter is straightforward. 

The analytical result contains logarithmic functions of rather complicated argu- 
ments that we schematize as 

a + PRi +- f R 2 + SR 1 R 2 NjRufy) 
° g a + (3R 1 -<yR 2 -5R 1 R 2 ° g N(R U -R 2 ) ' [ ' 

where a, /3, 7, S are polynomials and Ri and R 2 are square roots of polynomials of 
the radial variables; we note that the logarithms appearing in the result obtained using 
the expansion in Gegenbauer polynomials have a similar complicated form. In both 
cases, for continuing the analytical integration over the radial variable \q\ it is essential 
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to decompose these logarithmic functions [5] , by separating the various combinations of 
square roots Ri so that only functions of the form 

. a + bR 1 c + dR 2 , e + fR x R 2 

lo § — ro~ ' log — Td~ ' log — p 

a — oRi c — dR 2 e — } R\R 2 
appear. In fact, by observing that the product 

N(R 1 ,R 2 )N(-R 1 ,R 2 )N(R 1 , -R^)N(-R U -R 2 ,) 
is a polynomial squared and that 

N(Ri, R 2 )N(-Ri,R 2 ) = P 1 (c + dR 2 ) 2 , 
N(R ll R 2 )N(R 1 , -R 2 ) = P 2 (a + bR^ 2 , 
N(R 1 , R 2 )N(-R U -R 2 ) = P 3 (e + fRtR 2 ) 2 , 
(Pi are polynomials), we can rewrite eq.(8) as 



log 



so that 



N(R 1 ,R 2 ) _ 1 N(R 1 ,R 2 )N(-R 1 ,R 2 ) N(R U R^Nj-Ri, -R 2 ) 
N(R U -R 2 ) 2 °^ N(Ri, —R 2 )N(—Ri, —R 2 ) N(R 1: —R 2 )N(—Ri,R 2 ) ' 

. N(R 1 ,R 2 ) c + dR 2 e + fR 1 R 2 

log — — — = log — — + log 



N(R U -R 2 ) ^ c-dR 2 ' e- fR x R 2 ' (9) 
Applying the decomposition here exemplified to all the logarithmic functions, the 
analytical expression of the integral (4) takes the simpler form: 



/ 



da 
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[(q - fci) 2 + mj][(q - k 2 f + m 2 ] q^R(kl k 2 , (h - k 2 ) 2 ) 
Ci + R(k 2 l7 kl (h - k 2 ) 2 )R(q 2 , k 2 , -to 2 ) 



x 




d - R(kf 
C 2 + R(k 2 



C 2 - R(k 2 
C 3 + R(k 2 



C 3 - R{k( 



kKh-k^Riq^kl-m 2 ) 
kKh-hrW^i-m 2 ) 



kKh-k^Riq^kl-m 2 ) 
kl(h-k 2 ) 2 ' 



(10) 



kl (h - k 2 ) 2 ) 



where 



Ci =(k 2 - kl - (h - k 2 ) 2 )q 2 + kf + m\(k\ + k 2 - (h - k 2 f) 

- k 2 (kl + (h - k 2 ) 2 + 2m 2 2 ) , 
C 2 =(k 2 -k 2 - (h - k 2 ) 2 )q 2 + kj + m\(k\ + k 2 - (h - k 2 ) 2 ) 

-kKki + ih -k 2 ) 2 +2ml) , 
C 3 =(h ~ k 2 ) 2 -k\-k\. 



(11) 



The symmetry of eq.(10) under the exchange (fci,mi) <-> (£2,7712) is evident; we 
underline that eq.(10) has a form simpler than the result obtained using the expansion 
(2), which contains more complicated logarithms of the kind (8). 

Let us continue our analysis, considering this time the angular integral containing 
three denominators 

^ (12) 



[fo _ fcl )2 + m?] [(g _ k2 )2 + m 2] [( ? _ fc 3 )2 + m§] 

which is the hyperspherical angular part of the amplitude of the box graph shown in 
Fig. 2. Using again the expansion (2), we work out a triple series containing the angular 
integrals 

' Cl{q ■ k) C*(q ■ k 2 ) Ci(q ■ h) dSl q . (13) 



The analytical expressions of above integrals are not known; the known summation 
formulae for hyperspherical harmonics, which allows to work out eq.(6), are not useful 
in this case; the reason of this behaviour will be more clear in the following. Therefore 
the integral (12) cannot be calculated in analytical form using the method shown in 
ref. [1] . Thus, following the above described two-denominators example, we combine the 
three denominators of eq.(12) introducing two Feynman parameters, we perform the 
angular integration using the formula 



/ 



dn q 2** (14) 



[(q - k) 2 + m 2 ] 3 R 3 (q 2 ,k 2 ,-m 2 ) 



and the subsequent integrations over the Feynman parameters which are straightforward 
but rather cumbersome; the result is a sum of four logarithmic functions of asymmet- 
rical huge arguments which are decomposed with the same mechanism of eq.(8-9) and 
recombined into three logarithms of simpler arguments. We expect the result to have a 
structure similar to that of eq.(10). In fact one obtains 



/ 



dQ a it 2 

: X 



[(q-k 1 ) 2 + m 2 ] [(q - k 2 ) 2 + m 2 ] [(q - k 3 ) 2 + m 2 ] 



d + Rjq^kl-mpy/A ^ C 2 + R(q 2 , kj -m 2 ) VA 
1 ° g Ci - R(q 2 , k 2 , -m 2 )^A ° g C 2 - R(q 2 , k 2 , -m 2 )VA 

C 3 + R(q 2 ,k 2 ,- m 2 )VA 
" , C 3 -R(q 2 ,k 2 ,-m 2 )VA f 

where A is the polynomial of degree 2 in q 2 

A = a 2 (q 2 ) 2 + aiq 2 + a 
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with the coefficients 
a 2 =R 2 (a 1 , a 2 , a 3 ) 

«i =2wi(2ai6i - a\b 2 - a\b 3 + a 2 6 3 + a 3 b 2 - 02^2 - a 3 b 3 ) 

4 

+ 2(tti — b\)a\{a\ — a 2 — a 3 ) — -a\a 2 a 3 (+2 cyclic permutations) 

a =u\{a\ - 2a\b 2 - 2a\b 3 - 2b 2 b 3 + b\ + &§) 

+ 2mu 2 (-aia 2 + a\b\ + a\b 3 + a 2 b 2 + a 2 b 3 - 2a 3 b 3 - bi b 2 + b\b 3 + b 2 b 3 - b 3 ) 
+ 2ui{a\bi - a\a 2 b 2 - a\a 3 b 3 - a\b\b 2 - a\b\b 3 + 2aib 2 b 3 + a 2 b\ - a 2 b 2 b 3 
— a 3 b 2 b 3 + (X3&3) + a\b\ — 2a\a 2 b\b 2 (+2 cyclic permutations) 

where a t = (k 2 - k 3 ) 2 , a 2 = (ki - k 3 ) 2 , a 3 = {ki - /c 2 ) 2 , bi = kf and Ui = mj (we 
show a partial list of the expressions of «o an d ol\ with the indication that the complete 
expressions must be obtained summing up also the two cyclic permutations of the indices 
1,2 and 3 of the shown portions). 

The arguments of the logarithms contain 

Ci =f3 2 q 2 2 +[3iq 2 + A) 
f3 2 =ai - a 2 - a 3 

(3\ =u\(2a\ - a 2 - a 3 + 2b\ - b 2 - b 3 ) + u 2 (-a 2 - 61 + 63) + u 3 (-a 3 - bi + b 2 ) 
- 2ai6i - 2a 2 a 3 + a 2 b\ + a 2 b 2 + a 3 bi + a 3 b 3 

(3 G =u\(a\ -b 2 - b 3 ) + uiu 2 (-a 2 + bi + b 3 ) + uiu 3 (-a 3 + bi + b 2 ) 
+ u 2 bi(-a 2 + bi - b 3 ) + u 3 bi(-a 3 + 61 - b 2 ) + u 2 u 3 (-2bi) 
+ ui(2aibi - a 2 b 2 - a 3 b 3 - bib 2 - bib 3 + 2b 2 b 3 ) + aibf - a 2 b\b 2 - a 3 b 1 b 3 

and 

C 2 =Ci(with the cyclic permutation 1 — > 2, 2^3, 3^1) 
C 3 =Ci(with the cyclic permutation 1-^3, 2 — > 1, 3^2) . 
We note that: 

1) Due to the effect of decomposition (8-9), eq.(15) has a form such that the sym- 
metry under the permutations of the indices 1,2 and 3 is evident at sight. 

2) Comparing eq.(15) and eq.(10) we observe that the square root VK plays a role 
analogous to that of R(kf , k 2: (k\ — k 2 ) 2 ) in eq.(10), but it has a more complex structure, 
not reducible to one single square root of the kind R(x, y, z). This fact explains why the 
expansion in hyperspherical harmonics fails to obtain our result eq.(15): if the integrals 
of three Gegenbauer polynomials (13) were known in analytical form, we would be 
able to perform the angular integration of the product of the expansions of the three 
denominators of eq.(12); but the result would be the expansion of our result (15) in 
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Gegenbauer polynomials, surely very complicated, whose summation in analytical form 
to (15) would be almost hopeless. 

3) The angular integral of the box graph (15) can be used as a "building block" in 
order to work out more complex angular integrals. 

3. In this section we apply the result of the previous section to the analytical calculation 
of the contribution to electron g-2 of a typical "triple-cross" graph. 

We consider the scalar amplitude / of the not planar self-energy graph obtained 
from the vertex graph of Fig.l, disregarding the external photon interaction and the 
numerators of the Feynman integral. 

Following the section (2), we decompose the integral in radial and angular parts: 

l" dk\ r kl dkl r q 2 dq 2 2 2 2 2 

J ~J ~2~J %2TT) ^ 1 2 ,p) 

where J(k 2 , kl, q 2 ,p 2 ) is the angular integral 

7(1,2 7,2 2 2\ _ f dVlk 2 f dQ q 

J{k lt k 2t q,p)-J {p _ k2)2 J {q ^ k2) 2 [{p + q _ k2) 2 + 1] X 



X 



! d Jhl (16) 

J (h-p) 2 [(h-k 2 ) 2 + l] [fa-qF + l] { } 



and the electron mass has been set to 1. We underline that this graph is not planar 
and, from a topological point of view, it is made up of three box loops with one or 
two common sides; the momentum routing has been chosen in such a way that only 
one denominator be a linear combination of three momenta (this fact is unavoidable 
because of the not planarity of the graph). As we have shown in the previous section, 
the expansion of the denominators is not useful for evaluating this angular integral in 
analytical form, therefore we apply eq.(15) to the angular integration over ki, so we 
find: 



J (k 2 ,2 2, _ *f f dVL k2 r dQq ^ 

J^^q )- g2 J {p _ k2? ] (q-k 2 ) 2 [(p+q-k 2 ) 2 + l]£i 



'' logi^ 



(q-hrHp + q-hr + l}^ v/A ' 

the arguments of the logarithms contain the square roots R(kf,p 2 , —1), R(kf, kl, — 1), 
R(ki, q 2 , — 1), VA, whose only \/A depends on the angles among p, q and k 2 . We 
introduce hyperspherical angles for the remaining angular integrations 

f r 2 r 2 f 2n 

ditk 2 / dfl q = / 47rsin 9id6i / 27rsin (pid(pi / sin 4> 2 d(f) 2 
J Jo Jo Jo 
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where 9\ and are respectively the angle between p and k 2 and between q and 
and (j>2 is the azimuthal angle between p and q; the result of the integration over 02 is 
schematized as 

J(h 2 ,2 2n 4 r sin 2 gi r sin 2 0! ^ 

J^^J-tt ^ p2 + jb 2_ 2^11^1 cos (JiVo ? 2 + fc 2 2 -2| 9 ||/ C2 |cos0/ 0lX 

1 30 

X-=£Li 2 (G0 ( 17 ) 

VAo i=1 

where A is a polynomial of degree 2 in cos6*i and cos0i. A careful analysis of the 
structure of the square root \fKo and of the zeroes of the two denominators of eq.(17) 
containing cos 9\ and cos 0i shows that the angular integration over 9i is of elliptic kind, 
therefore the angular integral (16) cannot be reduced in closed analytical form. On the 
contrary we have worked out a integral representation. This not expected behaviour 
persists even using different momentum routings; it is probably due to the not planarity 
of considered graph. 

However, this fact does not stop our analytical calculation. We found that the 
total angular integral of the two-loop self-energy graph with two photon crossed is 
expressed by a elliptic integral involving logarithms obtained from eq.(8) if one uses a 
particular momentum routing; in this case an exchange between the radial integrations 
and the remaining angular integral avoids the ellipticity and allows one to complete the 
analytical calculation. 

Considering our three-loop result (17), one finds that Aq is a polynomial of degree 
two in q 2 , so that an exchange between the integrations over 9\ and q 2 allows one to 
continue the calculation. Work is in progress to extend this mechanism to the remaining 
radial integrations. 

All the algebraic manipulations were carried out through the symbolic manipulation 
program ASHMEDAI [7]. 
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Figure captions 

Fig.l: An example of triple-cross vertex graph. The cross indicates the external photon 

interaction. 

Fig. 2: The box graph. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9404203vl 



This figure "figl-2.png" is available in "png" format from: 
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